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Abstract 



On the basis of a simple exactly solvable model we discuss the possibilities for 
C^ ' state preparation and state control of atoms in a periodic optical potential. 

CJ ' In addition to the periodic potential a uniform force with an arbitrary time 

J>^' dependence is applied. The method is based on a formal expression for the full 

^^- evolution operator in the tight-binding limit. This allows us to describe the 

dynamics in terms of operator algebra, rather than in analytical expansions. 
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The energy eigenvalues of a quantum particle moving in a periodic potential form energy 

^ , bands (the Bloch bands) that are separated by bandgaps. The eigenstate within a band is 

jS^ , characterized by the quasimomentum, which determines the phase difference between two 

Ph. 

points separated by a period. An initially localized wavepacket typically propagates through 

r> I space, leading to unbounded motion. When an additional uniform force is applied, the Bloch 

bands break up into a ladder of equally spaced energy levels, which are called the Wannier- 

Stark ladder. In this case, a wavepacket of the particle extending over several periods can 

exhibit bounded oscillatory motion, termed Bloch oscillation, at a frequency determined by 



the level separation in the ladder. These early results of the quantum theory of electrons 
in solid crystals [1-4] have regained interest recently due to the advent of optical lattices 
for atoms. These lattices are formed when cold atoms are trapped in the periodic potential 
created by the superposition of a number of traveling light waves [5-8] . In contrast to the 
case of electrons in crystal lattices, these optical lattice fields have virtually no defects, they 
can be switched on and off at will, and dissipative effects can be largely controlled. The 
phenomenon of Bloch oscillations has first been observed for cesium atoms in optical lattices 
[9]. The uniform external force is mimicked by a linear variation of the frequency of one 
of the counterpropagating traveling waves, thereby creating an accelerated standing wave. 
By applying a modulation on the standing-wave position, Rabi oscillations between Bloch 
bands, as well as the level structure of the Wannier-Stark ladder has been observed for 
sodium atoms in an optical lattice [10]. Theoretical studies of transitions between ladders 
have also been presented [11]. Bloch oscillations have also been demonstrated for a light 
beam propagating in an array of waveguides, with a linear variation of the refractive index 
imposed by a temperature gradient [12]. 

When the applied uniform force is oscillating in time, the motion of a particle in a 
periodic potential is usually unbounded. However, it has been predicted that the motion re- 
mains bounded for specific values of the ratio of the modulation frequency and the strength 
of the force [13]. Similar effects of dynamical localization, including routes to chaos, have 
been studied experimentally for optical lattices, including both amplitude and phase mod- 
ulation of the uniform force [14]. Phase transitions have been predicted for atoms in two 
incompatible periodic optical potentials imposed by bichromatic standing light waves [15]. 

In the present paper we discuss the Wannier-Stark system with a time-dependent force, 
as a means of preparing the state of particles in a periodic potential. We derive an exact 
expression for the evolution operator of the particle, with an arbitrary time-dependent force. 
This allows one to apply the combination of delocalizing dynamics in the absence of the 
uniform force with the periodic dynamics induced by a uniform force for coherent control of 
the state of the particles. Exact solutions in the case of a constant uniform force have been 



obtained before by analytical techniques [17,18]. The operator method phenomena induced 
by an oscillating force to be described exactly in a unified scheme. Examples are dynamical 
localization and fractional Wannier-Stark ladders. 

The model is described in one dimension. However, this is no real restriction. Under the 
assumption of nearest-neighbor interaction, the corresponding 2D or 3D problem exactly 
factorizes into a product of ID solutions. 

II. MODEL SYSTEM 

A. Periodic potential 

The quantummechanical motion of atoms in a periodic optical potential V{x) with period 
a, is described by the Hamiltonian 

ffo = ^ + V-(.) . (1) 

We assume that the atoms are sufficiently cooled, so that only the lowest energy band is 
populated. The ground state in well n located at a; = na is indicated as \n). These states 
play the role of the basis of localized Wannier states. For simplicity we make the tight- 
binding limit, where only the ground levels in neighboring wells are coupled. When we 
choose the zero of energy at the ground level in a well, the Hamiltonian (1) projected on 
these ground levels is defined by 

Ho = -hQ{B+ + B_) , B±\n) = \n±l). (2) 

The raising and lowering operators i?+ and B^ are each other's Hermitian conjugate, and 
each one of them is unitary. The frequency fi measures the coupling between neighboring 
wells, due to tunneling through the barriers. We shall allow the coupling to depend on time. 
The eigenstates of Hq are directly found by diagonalizing the corresponding matrix. These 
states are the Bloch states |A;), with energy E{}z) = hQcos{ka). Their expansion in the 
Wannier states, and the inverse relations can be expressed as 



^ ^ inkai^\ |„\ / I Jl„ „—inka\ 



l^) = V^5^ 1^) ' l^^ = V2^i dke-^-''^\k) . (3) 

Obviously, the states \k) are periodic with period 27r/a, and the quasimomentum k can be 
chosen from the Brillouin zone [— 7r/a, ir/a]. The integration in (3) extends over this Brillouin 
zone. From the translation property {x\n) = {x + a\n + 1) of the Wannier wave functions 
it follows that the states (3) do indeed obey the Bloch condition {x + a\k) = exp{ika){x\a). 
When the states \n) are normalized as {n\m) = Snm, the Bloch states obey the continuous 
normalization relation {k\k') = 6{k — k'). 

B. Uniform force 

An additional uniform force is described by adding to the Hamiltonian the term 

i?. = — . (4) 

a 

where the (possibly time-dependent) force of size hA{t)/a is in the negative direction. On 
the basis of the Wannier states, this term is diagonal, and it is represented as 

Hi = hABo , Bo\n) = n\n) . (5) 

Hence the evolution of a particle under the influence of the total Hamiltonian 

H = Ho + Hi, (6) 

with Hq and Hi defined by eqs. (2) and (5), in terms of the operators B± and Bq. We shall 
also need expressions for the operators B± and Bq acting on a Bloch state. These can be 
found from the definition of the operators and the expansions (3). One easily finds that 

B±\k) = e^^'^'^IA;) , e-'^^»|A;) = \k - -) . (7) 

a 

In Bloch representation the operators have the significance B± = exp(=p?A;a), Bo = 
{i/a){d/dk), which is confirmed by the commutation rules (8). The Wannier states may be 
viewed as discrete position eigenstates, with Bq the corresponding position operator. The 
Bloch states play the role of momentum eigenstates, and the finite range of their eigenvalues 
within the Brillouin zone reflects the discreteness of the position eigenvalues. 



C. Operator algebra 

The basic operators B±. and Bq obey the commutation rules 

[Bo,B^] = ±B^, [B+,B_]=0. (8) 

In order to derive exact expressions for the evolution operator corresponding to the Hamil- 
tonian (6), we need several operator identities involving these operators Bq and i?±. The 
identities 

ei/3i?o5^g-i/3i?o = e^'^B± . (9) 

directly follow from the commutation rules (8), and they lead the transformation rules 

^i^Bo g^p (^ _ i-a{B+ + B_))e-'^^'-' = exp ( - i-a{e'^B+ + C'^B^)) (10) 

for arbitrary values of a and 13. We shall also need the equalities 

% t t 

exp{-aB±)Boexp{--aB±) = BqT 2^B± , (11) 

which are verified after differentiation with respect to a, while using the commutation rules 



III. OPERATOR DESCRIPTION OF EVOLUTION 

A. Evolution operator 

In this section we derive expressions for the evolution operator U{t, 0), which transforms 
an arbitrary initial state |^I/(0)) as |^(t)) = U{t, 0)|^(0)). The results are valid for any time- 
dependence of the uniform force and the coupling between neighboring wells, as specified 
by A(t) and Q{t). A time-dependent coupling represents the case that the intensity of the 
lattice beams is varied. We express the evolution operator in the factorized form 

[/(t,0) = f/i(t,0)f/o(t,0), (12) 



where Ui{t,0) = exp[—i(f){t)Bo] gives the evolution corresponding to the Hamiltonian Hi 
alone, in terms of the phase shift 

0(t) = /* dt' A(t') . (13) 

From the evolution equation for U with the Hamiltonian (6), while using the transformation 
(9) we find the evolution equation 

dUo ^ _^M(e^<AW5^ + e-'f^'^B)Uo{t) . (14) 

Since this equation only contains the commuting operators 5+ and i?_, it can easily be 
integrated. In fact, the solution is given by eq. (10) with the time- dependent values of the 
real parameters a and /3 defined by the relations 

a(t)e*'3(*) = /* rft'fi(t')e''^(*') . (15) 

JO 

Combining this solution with the definition of f/i, leads to a closed expression for the evolu- 
tion operator f/(t, 0) for an arbitrary time dependence of the uniform force, in terms of the 
parameters a, /5 and 0, defined in (13) and (15). The result is f/(t, 0) = i?(a,/9, 0), with R 
defined by 

i?(a, /3, 0) = e^(^-'^)^« exp ( - i\[B^ + 5_))e-^^^« . (16) 

This defines the unitary operator i? as a function of the three parameters a, /3 and 0. The 
result is valid for an arbitrary time dependence of the force and the coupling, described by 
A(t) and r2(t). The characteristics of the evolution of an arbitrary initial state is determined 
by the properties of the operators i? as a function of a, /9 and 0. Mathematically, these 
operators form a three-parameter group, which is generated by the three operators fi-t and 

5o. 

On the basis of the Wannier states, the contribution of the operator Bq in (16) is trivial, 
whereas the effect of the exponent containing B± can be evaluated by first expanding a Wan- 
nier state in Bloch states, for which the action of this exponent is simple. Then reexpressing 
the Bloch states in Wannier states, we find 
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exp ( - i]-a{B+ + B_)) \m) = Y. r"+'" J„_^(a) \n) , (17) 

n 

where we used the defining expansion exp(i^sin0) = Z]nexp(m0) J„(,^) of the ordinary 
Bessel functions. Hence the matrix elements of the operator (16) between Wannier states 
are 

{n\R{a,P,(P)\m) = (^e-^'^)-"+'"e-^"^J„_„(a) . (18) 

For the evolution operator (16) in Bloch representation we can just use the form of the 
operators B± and Bq, as given in Sec. II B. This leads to the result 

i?(a,/5,0)|A;) = e-'"™^(*^'^-^)|A;-0/a) . (19) 

This shows that the quasimomentum as a function of time varies as k{t) = k{0) — (f){t)/a, 
with 0(t) given in (13). The parameter determines the shift of the quasimomentum during 
the evolution. The expressions (18) and (19) clarify the significance of the three parameters 
a, (3 and that specify the evolution operator. 

B. Heisenberg picture 

The transport properties of any initial state is conveniently described by the evolution 
of the operators in the Heisenberg picture. Since any evolution operator can be written 
in the form of /?(«, /3, 0) for the appropriate values of the parameters, we can view R^BR 
as the Heisenberg operator corresponding to any operator B. The Heisenberg operators 
corresponding to B^ can be expressed as 

R^{a, (3, (p)B±R(a, (3, 0) = e^''^B± , (20) 

which is directly shown by using eq. (9). Since B± = exp{^ika) in Bloch representation, 
this confirms the significance of as the shift of the value of the quasimomentum. 

After using the transformation property (11), one finds the Heisenberg operator corre- 
sponding to the position operator Bq as 

7 



R\a,(3,<l))BoR{a,(3,<l)) = Bo + -^^i^-'^B. - e'^B+) . (21) 

This implies that the expectation value of the position after evolution is determined by 

{n) = {Bo) + -{e-'^^{B^)-e^^{B^)), (22) 

where the averages in the r.h.s. should be taken with respect to the inital state. Hence no 
displacement of a wavepacket can occur whenever (5+) = (-B_)* = 0. This is true whenever 
the initial state is diagonal in the Wannier states \n). Conversely, average motion of a 
wavepacket can only occur in the presence of initial phase coherence between neigboring 
Wannier states. The width of a wavepacket is determined by the expectation value of the 
square of the Heisenberg position operator (21). This gives the expression 

{n^) = {B^,) + ^{2-e-'^^{Bl)-e'^^{Bl)) 

+ y (e-^^(5ofi- + B^Bo) - e'^{BoB+ + B^Bo)) , (23) 



IV. LOCALIZED INITIAL STATES 

A. Arbitrary wavepackets 

A fairly localized initial state |^&(0)) = I]nC„|n) with a reasonably well-defined quasimo- 
mentum can be modelled by assuming that neigboring states have a fixed phase difference 
9, so that 

c„Cn+i = |c„c„-|_i|e . (24) 

Then the quasimomentum is initially centered around the value /cq = 0/a- For simplicity, we 
assume moreover that the distribution over Wannier states is even in n, so that |c„| = |c_„|. 
The initial average position of the particle is located at n = 0. In order to evaluate the 
time-dependent average position and spreading of the packet, we can apply eqs. (22) and 
(23). The symmetry of the distribution implies that (Bq) = 0, while (B^) = a^ is the initial 
variance of the position. When we introduce the quantities 



Y^ |C„+1C„| =bi , Y^ \Cn+2Cn\ = h , (25) 

n n 

we obtain the simple identities 

(5+) = he-'' , {BD = he-''' , (5o5+) = -{B+Bo) = h,e-'' . (26) 

The last identity is proven by using that the quantity /2n+i = |c„+ic„| is even in its index 
(which takes only odd values). Therefore, J2iU'i = 0, which is equivalent to the statement 
that 2{B^Bq) + (-B+) = 0. The other expectation values occurring in eqs. (22) and (23) 
are found by taking the complex conjugates of the identities (26). This leads to the simple 
exact results 

{n) = abi sin(/3 - 9) , (n^) = a^ + —[l-b2 cos 2(/5 - 9)) , (27) 

so that the variance of the position is found as 

^2 ^ ^^2^ _ ^^^2 ^ ^2 ^ ^(l _ 52 _ (^^ _ ^2) ^^g2(/3 - 9)) . (28) 

Notice that the parameters 61 and 62 are real numbers between and 1. In the limit of 
a wide initial wavepacket, determined by coefficients c„ whose absolute values vary slowly 
with n, the parameters 61 and 62 both will approach 1, and the width a will not vary during 
the evolution. In the opposite special case that the initial state is the single Wannier state 
|0), one finds that 61 = 62 = 0, so that the width a = a/\^. 

In the special case that the particle is initially localized in the single Wannier state at 
X = 0, so that I^E'(O)) = |0), the parameters 61, 62 and (Tq vanish, so that 

(n) = , a^ = (n^) = 0^2 . (29) 

This shows that the average position of the wavepacket does not change, and that its width 
is determined by the parameter a alone. This is in line with the fact that the population 
distribution over the Wannier states after the evolution is pn = |(n|i?|0)p = J^(a), as 
follows from Eq. (18). Hence the (time-dependent) value of a determines the spreading of 
an initially localized particle. 



B. Gaussian wavepacket 

When the initial distribution over the sites is Gaussian with a large width, we can 
evaluate the full wavepacket after evolution. Suppose that the initial state is specified by 
the coefficients 

e^-'eM-^). (30) 



which obey the condition (24). This state is properly normalized provided that a^ ^ 1. 
When the evolution operator is expressed as in (16), the time-dependent state is expanded 
as |\I'(t)) = -R|^(0)) = Z]n/nexp[m(^ — 0)]|n). Summation expressions for the coefficients 
/„ are directly obtained by using the expression (18) of R in Wannier representation. We 
use similar techniques as applied in Ref. [16] in the context of the diffraction of a Gaussian 
momentum distribution of atoms by a standing light wave. The technique is based on 
differentation of the expression for /„ with respect to n, while using the property a[J„+i(Q;) + 
J„_i(a)] = 2nJn{<y) of Bessel functions. When the width is sufficiently large, so that the 
difference /„+i — /„ can be approximated by the derivative, this leads to the differential 
equation 



2 4/ n _, /_ „;„/,3 n\ „A x I ,•„ „//3 n\ J^ 



2af,^^ [asinil3-e)-n]fn + iacos{l3-e)^ . (31) 

an ^ ' an 

By solving this equation, we arrive at the closed expression 

/„ = -: exp ( — -^ : — — -— j , 32) 

A/ ^ 2(Jq — la cos(p — 0) ' 

with the normalization constant determined by 

^ = ^(2^0 + 1^^ -) ■ (33) 

We find that the distribution is Gaussian at all times, with a time-varying average 
position and variance. These are given by the expressions 

2 

(n) = a sin(/3 - 9) , a^ = a^ + ^(l + cos 2(/3 - 9)) . (34) 
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These results are in accordance with eqs. (27) and (28), as one checks by using the ap- 
proximate expressions bi = exp(— /^/8o"q) ^ 1 — P/ScTq, while neglecting terms of order 
(l/(To)^ and higher. The width of the packet never gets smaller than its initial value. The 
phase difference between neighboring sites is mainly determined hj 6 — (f). This shows that 
a phase difference can be created or modified in a controlled way, simply by imposing a 
time-dependent force that gives rise to the right value of 0. Notice that in these expressions 
(34), 6 and f3 enter in an equivalent fashion. The position and the width of the Gaussian 
distribution can be controlled at will by adapting the force to the desired value of f3. 

We recall that the results of this chapter are valid for an arbitray time-dependent force 
A(t), which determines the time-dependent values of the parameters a, (3 and as specified 
in eqs. (13) and (15). In the subsequent sections, we specialize these expressions for constant 
or oscillating values of the uniform force. 

V. CONSTANT UNIFORM FORCE AND BLOCK OSCILLATIONS 

A. Wannier-Stark ladder of states 

The case of a constant force is the standard situation where Bloch oscillations occur. 
When A and Q are constant, the Hamiltonian is time-independent, and then it is convenient 
to introduce the normalized eigenstates \ipm) of H. When we expand these eigenstates in 
the Wannier states as \^m) = Y.n l^)'^!™'-', the eigenvalue relation H\tjjm) = -^ml'^m) with 
Em = hujm leads to the recurrence relations for the coefficients 

l^{ct\ + ct\) + ^nc^r^ = ^rr^C ■ (35) 

We introduce the generating function 

^n.(A:) = ./f Ec^^e-'=^ (36) 

V 27r „ 

which is normalized for integration over the ffist Brillouin zone. In fact, from the expression 
(3) of the Bloch state, one notices that the generating function Zm{k) = (klipm) is equal to 

11 



the Bloch representation of the eigenstate \ipm)- The relations (35) are found to be equivalent 
to the differential equation 

VLcos{ka)Zm{k) - ——Zm{k) = UmZm{k) , (37) 

with the obvious normalized solution 

Zm{k) = W7— exp ( — [r2sin(A;a) — akujm]) . (38) 

Since the functions Zm{k) as defined by (36) are periodic in k with period 27r/a, the same 
must be true for the expressions (38). Hence, the frequency eigenvalues must be an integer 
multiple of A, so that we can choose Um = mA, with integer m. For these values of the 
eigenfrequencies, the coefficients c^™-* follow from the Fourier expansion of Zm, with the 
result 

Ci™) = {n\iJm) = Jn.-n{n/A) . (39) 

We find that the total Hamiltonian H has the same eigenvalues as Hi. Apparently, the 
energy shifts due to the coupling between the Wannier states as expressed by Hq cancel 
each other. Since the energy eigenvalues are integer multiples of A, each solution of the 
Schrodinger equation is periodic in time with period 2tt/A, and the same is true for the 
evolution operator U{t) given in eq. (16). This also implies that an initial localized state 
remains localized at all times, due to the addition of the uniform external force. The eigen- 
states \ipm) are the Wannier-Stark ladder of states [10]. They form a discrete orthonormal 
basis of the first energy band, and they are intermediate between the Wannier and the Bloch 
basis of states. 

B. Oscillations of localized states 

The definitions (13) and (15) show that 

a = (21]/A) sin(At/2) , p = At/2 , = At . (40) 

12 



In the Wannier representation, the matrix elements of U are found from (16) as 

(n|f/(t,0)|m) = r"+-e-'^*("+™)/V„_„(^sin^) , (41) 

which represents the transition amphtude from an initial state \m) to the final state \n). 
For the initial Wannier state |^I/(0)) = |0), the time-dependent state is |^E'(t) >= J2n fn{t)\n) 
with 

/„(t)=r"e-^^*"/V„(^sin^). (42) 

This is in accordance with Eq. (50) of ref. [17], which has been obtained by a rather 
elaborate analytical method, rather than an algebraic one. Equation (29) shows that the 
time-dependent average position (n) of the wavepacket remains zero at all times, whereas 
the mean-square displacement a = |a|/V2 displays a breathing behavior, and returns to 
zero after the Bloch period 2tt/A. Moreover, according to eq. (42), the phase difference 
between neighboring sites varies continuously with time. 

This is already quite different when only two Wannier states are populated initially. 
Consider the initial state 

\nO)) = ^m + e^'\l)). (43) 

Then the average position can be evaluated from eq. (22), for the values of a and f3 as given 
in (40). The result is 

(n) = ^ + ^(cos^-cos(At-^)), (44) 

which shows that the packet displays a harmonically oscillating behavior. The amplitude 
of the oscillation is governed by the ratio Q/A, which is half the maximum amplitude for 
Bloch oscillations of a wavepacket with a large width (see Sec. V C). This amplitude must be 
appreciable in order that interband coupling induced by the uniform force remains negligible, 
as we have assumed throughout this paper. The distribution pn = |/np after half a Bloch 
period, both for the initial single Wannier state and for the inital state (43) is illustrated in 

13 



Fig. 1. This demonstrates that a strong displacement can already be induced by evolution of 
a superposition state of just two neighboring Wannier states, with a specific phase difference. 
This displacement arises from the interference between the transition amplitudes from the 
two initial states to the same final state \n). 

C. Bloch oscillations and breathing of a Gaussian wavepacket 

The evolution of a Gaussian wavepacket as discussed in Sec. IV B is specialized to the 
present case of a constant force after substituting the expressions (40) in eqs. (32)- (34). We 
find for the average position (n) the identity 

(n(t)) = —[cos e - cos(^ - At)] . (45) 

This demonstrates that the wavepacket oscillates harmonically in position with frequency A, 
and with amplitude Q/A in units of the lattice distance a. The velocity of the wavepacket 
is found from the time derivative of (45), with the result 

v{t) = -ansm{e - At) . (46) 

It is noteworthy that this expression (46) coincides exactly with the expression for the 
group velocity dE/hdk, with the derivative evaluated at the time-dependent value of the 
quasimomentum {6— At) /a, with E = hQcos{ka) the dispersion relation between energy and 
quasimomentum in the absence of the uniform force, as given in Sec. II A. Apparently, the 
expression for the group velocity retains its validity also in the presence of the uniform force. 
Of course, the concept of Bloch oscillations of the wavepacket as a whole has significance 
only when the amplitude Q/A of the oscillation is large compared with the width a of the 
packet, which in turn must extend over many lattice sites. 

The time-dependent width a of the Gaussian packet is found from eq. (34) in the form 

(T^ = al + j-2X^(l - cos At) (l + cos(At - 2^)) . (47) 
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Hence the variance of the position deviates from its initial value by an oscillating term. 
The amplitude of this oscillation is governed by the ratio (r2/2Ao"o)^- The initial width is 
restored whenever one of the terms in brackets vanish. This happens twice during every 
Bloch period, except when 9 = 7r/2, when these two instants coincide. This combined 
breathing and oscillating behavior is illustrated in Figs. 2 and 3, for various values of the 
relatice phase 6. Notice that the oscillation is always harmonic with the Bloch frequency 
A. This is due to the simple form of the dispersion relation for the case of nearest- neighbor 
interaction. The time dependence of the variance is a superposition of terms with frequencies 
A and 2A. 

D. Zero external force 

In the absence of the external force, we can take the limit A ^ in the results of the 
previous subsections. In particular, this gives = /3 = 0, a{t) = Qt. Then the evolution of 
an initial Wannier state |^(0)) = |0) is given by 

|^(t)) = R\ij{0)) = ^rV„(fit)|n) , (48) 

n 

which shows that the free spreading of an initial Wannier state after a time t gives Wannier 
populations equal to j9„ = | J„(r2t)p [19]. The mean-square displacement increases hnearly 
in time, as o" = nt/^/2. This shows that the spreading is unbounded in the absence of an 
external force. The self-propagator po{t) decays to zero for large times. The phase difference 
between neighboring sites is ±7r/2 at all times. For only two coupled wells, the coupling 
would give rise to Rabi oscillations with frequency Q. Equation (48) can be viewed as the 
generalization to the case of an infinite chain of wells. 

For a Gaussian wavepacket with initial width ctq and initial quasimomentum determined 
by 6, expressions (45) and (47) take the form 

(n(t)) = -fit sin ^ ^ a"^ = al + — 2-(l + cos2^) . (49) 
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As one would expect in the absence of a uniform force, the group velocity takes the constant 
value V = —aQ sin 6, which leads to unbounded motion of the packet (except for ^ = 
or ±7r). Usually, the width increases indefinitely during he propagation. However, for the 
special values 6 = ±7r/2 the width is constant, and the packets propagates as a solitary wave. 
Notice that such a phase difference between neighboring Wannier states arises spontaneously 
when a single Wannier state spreads in the absence of a uniform force. 

VI. OSCILLATING FORCE 

Other situations of practical interest arise when the uniform force has an oscillating 
component. Examples are the coupling between the states in the Wannier-Stark ladder 
[10], and dynamical localization for special values of the amplitude-frequency ratio of the 
oscillation [13,14]. The situation of an oscillating force is also decribed by the operator 
description of Sec. Ill A. We give some results below. 

A. AC force only 

The situation of a harmonically oscillating uniform force can be expressed as 

A (t) = (5 cos (cut) , (50) 

so that = {6/uj)sm{ujt). Then according to (15) the parameters a and P are specified by 
the equalities 

ae'^ = ntJoi-) + nJ2 ^(-) — fe*""* - l) , (51) 

where we used the expansion defining the ordinary Bessel functions, given in section III A. 
The first term in (51) increases linearly with time, whereas the summation is bounded, 
and periodic in time with period T = 1t\ juj. The behavior of a and /3 as defined by (51) 
is quite complicated in general. However, for large times the value of a, and thereby the 
spreading of an initial Wannier state, is the same as in the absence of the uniform force, 
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with Q replaced by the reduced effective couphng QJq{S/uj). After one period T, the values 
of the parameters become simple, and we find l3 = (p = 0, a = QTJq{6/uj). The evolution 
operator U{T) during one period T is simply given by the operator R defined in (16), at 
these values of the parameters. The eigenstates of the evolution operator R = U{T) are 
simply the Bloch states \k). The eigenvalues can be expressed as the exp{—iS{k)T/h), with 

£{k) = hnJoi-) (52) 

CO 

the corresponding values of the quasienergy, which are strictly speaking only defined modulo 
hu. The quasienergy bandwidth is reduced by the factor Jq{6/uj), compared with the energy 
bandwidth in the absence of the uniform force. 

When the ratio 6/uj of the amplitude and the frequency of the oscillating force coincides 
with a zero of the Bessel function Jq, no unbounded spreading occurs, and an initially local- 
ized state remains localized at all times, with a periodically varying mean-square displace- 
ment. The quasienergy bandwidth is reduced to zero in this case. This effect of dynamical 
localization has been discussed before for electrons in crystals [13]. A related effect of an ef- 
fective switch-off of atom-field coupling occurs for a two-level atom in a frequency-modulated 
field when the ratio of the amplitude-frequency ratio of the modulation equals a zero of the 
Bessel function Jq. This effect, which leads to population trapping in a two-level atom, has 
recently been discussed by Agarwal and Harshawardhan [20]. 

B. AC and DC force 

A constant uniform force creates Wannier-Stark states with equidistant energy values. 
An additional oscillating force can induce transitions between these states. Therefore, we 
consider the force specified by 

A{t) = Aq + 5 cos{ujt) . (53) 

Then the value of the parameters 0, a and (3 are 



17 



0(t) = Aot + {6/uj) sinM) , ae'^ = 1]^ J„(i)— -A ^MAo+n.)* _ ^\ ^^^^ 

„ u ?(Ao + nuj) ^ ^ 

In general, each term in the summation is bounded and periodic, but the different periods 
can be incompatible. Moreover, whenever ^Q + nu = 0, the corresponding summand attains 
the unbounded form QtJn{S/uj). At such a resonant value of Aq, the spreading of an initially 
localized state becomes unbounded, and the particle becomes delocalized. This delocaliza- 
tion is suppressed again when the ratio 6/uj is equal to a zero of the corresponding Bessel 
function J„. This is a simplified version of the phenomenon of fractional Wannier-Stark 
ladders, which has recently been observed and discussed [21,22]. 

The quasienergy values are again determined by the eigenstates of the evolution operator 
f/(T) for one period of the oscillating force. This operator is equal to the general operator 
R defined in (16), with the parameters 

a = 2f]sin(AoT/2) J] U-)-^ , /3(T) = AoT/2 , 0(T) = AqT . (55) 

These expressions are correct whenever Aq + nuj is nonzero for all values of n. Since these 
values of the parameters can be directly mapped onto the values (40) specifying the evolution 
with a constant uniform force, also the eigenvectors and corresponding quasienergies are 
immediately found. The eigenvectors of R can be expressed as \^m) = Y,n l^)cl™\ with the 
expansion coefficients c^™-* = Jm-n{C)- Here the argument ( of the Bessel functions must be 
chosen as the sum 

c=nj:u-)-^-^, (56) 

„ uj Aq + nuj 

which replaces the simple argument Q/A in eq. (39). The eigenvalues of i? = U{T) are 
exp{—iSmT/h), with the discrete quasienergy values Sm = hm Aq (modulo huj). 

In the resonant case that Aq+Hquj = for some integer no, one summand in the expression 
for a and (3 is modified, as indicated above. When T = t, only this modified summand is 
nonzero, and the evolution operator U{T) = R for one time period is characterized by the 
values 



a = nrJn^^ , /5 = , = -27rno . (57) 

The eigenvectors of R are the Bloch states \k), and the corresponding quasienergy values 



are 



S{k) = nQJn,{-)- (58) 



VII. DISCUSSION AND CONCLUSIONS 

We have analyzed the Wannier-Stark system, which is characterized by the Hamiltonian 
(6), in terms of the operators B± and Bq. The present interest in this model arises from 
the dynamics of atoms in a periodic optical potential, with an additionally applied uniform 
external force. We adopted the tight-binding limit, which implied nearest-neighbor inter- 
action only. This gives rise to an explicit simple dispersion relation between energy and 
quasimomentum, which makes the model exactly solvable. From the commutation proper- 
ties of the basic operators we obtain eq. (16) for the evolution operator for an arbitrary 
time dependence of the uniform force, where the three parameters are defined in (13) and 
(15). As shown in Sees. IIIB and IV, the parameter determines the shift in the value 
of the quasimomentum, whereas a and (3 determine the evolution of the average position 
and the width of a wavepacket. A particle starting in a single Wannier state has a uniform 
distribution over the quasimomentum, and cannot change its average position, wheras the 
width of its wavepacket is simply measured by a. On the other hand, even when only two 
neighboring states are populated initially, the wavepacket can display an appreciable mo- 
tion. In Sec. IV B it is demonstrated that an initially Gaussian packet remains Gaussian 
at all times. This remains true when the initial state has a non-zero expectation value of 
the quasimomentum, which is described as an initial phase difference between neighboring 
Wannier states. 

These results, which are valid for a uniform force with an arbitrary time dependence, 
unify and extend earlier results obtained for a constant or an oscillating uniform force. A 
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constant force induces Bloch oscillations of a wavepacket, and we obtain a simple expression 
for the amplitude of the oscillation and for the time dependence of the width of a wavepacket. 
For an oscillating force, the operator method shows that the quasienergy bands can be 
evaluated directly in terms of the value of the parameter a after one oscillation period. 
This produces an exactly solvable model for dynamical localization and fractional Wannier- 
Stark ladders. In general, by selecting a proper time-dependence of the force or of the 
coupling between wells, thereby realize the desired values of the parameters a, {3 and 0, 
we can coherently control the width and the position of a wavepacket, as well as the phase 
difference between neighboring sites. 
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FIGURES 
FIG. 1. Upper part: plot of the breathing population distribution for an initial Wannier state 

|0). Lower part: plot of the oscillating population distribution, for two initial superposition of 

Wannier states |0) and |1), and two different values of the relative phase 9. Both plots are evaluated 

for r^/A = 6. Shaded distributions hold after half a Bloch period t = vr/A. 

FIG. 2. Periodic behavior of the width and the average position of a Gaussian wavepacket for 
various initial values of the phase difference 6 between neighboring states. Initial value of the width 
is o-Q = 4, and l^/A = 50. 

FIG. 3. Bloch oscillation and corresponding breathing behavior of a Gaussian wavepacket in a 
constant uniform force. Values of (Tq, and and A as in Fig. 2. Upper part: = 0. Lower part: 
e = 7r/2. 
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